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SUMMARY 


A  digital  computer  method  of  determining  the  optimum  design  para¬ 
meters  of  lightly  loaded  shrouded  propellers  is  presented.  For  such  a 
case  the  geometry  of  the  vortex  wake  may  be  ideally  described  as  a  heli- 
cal  sheet  of  equal-pitch  vortex  filaments  arising  from  each  blade  trail¬ 
ing  edge  and  a  right  circular  boundary  sheet  composed  of  equal-pitch 
helical  filaments  arising  at  the  shroud  trailing  edge .  The  optimum  con¬ 
dition  for  light  loading  requires  that  this  wake  vortex  pattern  move  as 
a  rigid  body.  This  pattern  of  distributed  vorticity  is  arbitrarily  re¬ 
placed  with  a  finite  n;miber  of  finite,  but  unknown  strength,  vortex  fila¬ 
ments.  Using  the  Biot-Savart  relation,  an  expression  for  the  velocity 
components  associated  with  each  filament  is  obtained.  Applying  the 
boundary  conditions  to  these  sets  of  equations  results  in  a  set  of  si¬ 
multaneous  equations  in  terms  of  the  unknown  filament  strengths .  The 
blade  bound  vortex  distribution  obtained  from  the  simultaneous  solution 
is  compared  with  the  measurements  previously  obtained  from  potential  tank 
models  of  a  two-  and  four-bladed  shrouded  propeller.  It  is  concluded 
that  this  digital  computer  method  yields  results  that  are  adequate  for 
design  purposes,  and  it  is  recommended  that  this  method  be  used  instead  of 
the  potential  tank  method  for  obtaining  the  necessary  charts  and  tables 
covering  the  range  of  design  parameters  for  shrouded  propellers. 
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COWCUJSIONS 


The  digital  computer  method  of  determining  the  optimum  blade  bound 
vortex  strength  distribution  of  shrouded  propellers  presented  herein 
yields  results  very  similar  to  those  obtained  from  the  potential  tank 
models.  Of  the  two  approaches  to  the  problem,  it  is  concluded  that  the 
digital  computer  method  is  the  more  accurate  and  efficient.  The  use  of  a 
Fourier  series  representation  of  the  blade  bound  vortex  distribution  was 
proven  to  be  advantageous,  primarily  as  an  aid  in  presenting  the  results 
for  design  purposes. 

Although  the  method  was  developed  for  the  lightly  loaded  case ,  it 
is  believed  that  the  approach  with  certain  modifications  to  the  boundary 
conditions  may  be  employed  to  deteimine  the  optimum  blade  bound  vortex 
strength  for  heavily  loaded  shrouded  propellers. 
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RECOMMENDATIONS 


It  is  recommended  that  the  digital  computer  method  of  determining 
the  optimum  blade  bound  vortex  strength  distributions  of  shrouded  propel¬ 
lers  be  used  Instead  of  the  potential  tank  method  in  any  extension  of  the 
present  work.  In  addition^  the  use  of  a  Fourier  series  representation  of 
the  bound  vortex  strength  will  reduce  the  number  of  charts  and  tables  re¬ 
quired  to  cover  the  range  of  design  variables. 

Inasmuch  as  the  present  method  was  developed  for  lightly  loaded, 
shrouded  propellers,  it  is  believed  that  a  heavily  loaded  theory  would  be 
more  advantageous .  It  Is  recommended  that  any  extension  of  the  present 
work  be  for  the  heavily  loaded  case . 


5 


INTRODUCTION 


The  mathematical  desci^iptlon  of  fluid  flow  fields  associated  with 
lifting  surfaces  has  often  been  made  less  difficult  by  the  Introduction 
of  the  bound  vortex  concept  and  its  accompanying  pattern  of  vortex  sheets . 
This  has  been  particularly  true  for  those'  cases  of  three-dimensional  wings 
and  for  propellers  in  which  the  geometry  of  the  vortex  patterns  could  be 
ascertained  beforehand.  If  the  geometry  is  known  and  the  vortex  filament 
and  sheet  strengths  are  unknown^  then  the  solution  is,  in  theory, 
straightforward,  and  two  classical  avenues  may  be  followed.  One  method 
involves  the  obtaining  of  a  solution  to  Laplace's  equation  in  terms  of  the 
velocity  potential  which  satisfies  the  boundary  conditions  of  the  particu¬ 
lar  problem.  By  further  consideration  of  the  direct  analogy  that  exists 
between  the  electrical  and  velocity  potentials,  potential  tank  models  may 
be  used  to  obtain  solutions  for  those  cases  which  are  not  mathematically 
tractable.  In  addition,  of  course,  a  solution  may  be  reached  through  nu¬ 
merical  procedures  and  approximations.  A  second  method  makes  use  of  the 
Biot-Savart  equation  which  is  a  differential  relation  between  the  incre¬ 
mental  velocity  at  a  point  associated  with  an  incremental  length  of  vortex 
filament  and  its  strength  and  geometry.  The  integration  of  this  expression 
usually  Involves  numerical  procedures  and  approximations  and  also  the 
satisfaction  of  certain  boundary  conditions . 

In  recent  years,  a  number  of  prototype  and  test-bed  aircraft  hav¬ 
ing  VTOL/STOL  capabilities  have  been  constructed  that  have  used  the 
shrouded  or  ducted  propeller  for  propulsion  and  sustentation.  The  results 
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of  the  tests  of  these  aircraft  have  indicated  a  need  for  a  better  under¬ 
standing  of  the  effects  of  the  shroud  or  duct  on  propeller  performance. 
Several  theoretical  investigations  have  been  directed  toward  determining 
the  necessary  optimum  design  parameters  of  propellers  operating  in  ducts. 
One  analysis^  considered  a  lightly  loaded  propeller  having  a  minimum  ki¬ 
netic  energy  loss  operating  in  a  circular  duct  of  constant  diameter  and 
of  infinite  or  semi-infinite  length.  A  series-type  solution  of  Laplace's 
equation  for  the  velocity  potential  with  the  necessary  boundary  condi¬ 
tions  was  obtained.  The  coefficients  appearing  in  the  solution  were  de¬ 
termined  by  a  standard  method  of  successive  approximations.  The  bound 
vortex  strength  distribution  was  determined,  and  the  increase  in  this 
quantity  and  the  thrust  increment  over  that  of  a  free  propeller  were  pre¬ 
sented  in  tables  for  a  propeller  operating  in  a  duct  with  zero-tip  clear¬ 
ance.  The  tables  presented  the  values  for  three-,  four-,  five-,  and  six- 
bladed  shrouded  propellers  operating  at  advance  ratios  of  one-sixth  to 
two -thirds . 

The  potential  tank  technique  of  determining  the  optimum  design 

2 

parameters  of  shrouded  propellers  has  been  developed  and  presented  by 
the  author  of  this  report.  In  developing  the  technique,  potential  models 
were  constructed  of  the  ultimate  wake  of  a  two-  and  four-bladed  shrouded 

^A.  J.  Tachmindji,  The  Potential  Problem  of  the  Optimum  Propeller 
with  Finite  Number  of  Blades  Operating  in  a  Cylindrical  Duct,  David  Taylor 
Model  Basin,  Report  No.  1228,  July,  1958* 

^ .  B .  Gray ,  An  Investigation  of  an  Approach  to  the  Problem  of 
Determining  the  Optimum  Design  of  Shrouded  Propellers,  U.  S.  Army  Trans¬ 
portation  Research  Command,  Report  TREC  60-44,  May,  I96O. 
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propeller  operating  at  an  advance  ratio  of  1.556  and  the  measurements  were 
made  on^  and  the  data  presented  for,  these  models.  It  was  concluded  that 
the  approach  yielded  answers  suitable  for  design  pruposes  but  that  the 
model  construction  was  difficult  and  time  consimiing.  A  second  method  of 
determining  the  necessary  design  parameters  involving  the  digital  computer 
integration  of  the  Biot-Savart  relation  was  also  presented  in  that  report. 
Although  the  results  of  the  digital  computer  analysis  were  inconclusive  due 
to  limitations  of  time  and  funds,  it  was  recommended  that  the  computer 
approach  be  pursued  to  its  conclusions.  It  is  the  results  of  the  exten¬ 
sion  of  this  digital  computer  analysis  that  are  presented  in  this  report. 
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ANALYSIS  OF  PROBLEM 


For  the  case  of  the  lightly  loaded  shrouded  propeller,  the  geometry 

of  the  ultimate  vortex  wake  may  be  easily  described,  and  hence  the  well- 

3 

established  vortex  theories  may  be  applied  in  a  straightforward  manner. 
Here,  the  term  "lightly  loaded"  should  be  interpreted  as  yielding  in¬ 
duced  velocities  that  are  negligibly  small  in  comparison  with  the  free- 
stream  velocity  and  the  local  blade  rotational  velocities.  In  addition, 
it  is  implied  that  the  wake  has  reached  its  ultimate  configuration  at  the 
shroud  trailing  edge  and  that  a  vortex  filament  shed  at  a  nondimens ional 
blade  station  remains  at  that  nondimens ional  radius  with  respect  to  the 

local  wake  radius.  The  resulting  vortex  pattern  has  been  discussed  in 

2 

Appendix  I  of  Gray's  Report  and  its  geometry  is  assumed  to  be  speci¬ 
fied  by  the  free-stream  velocity  and  the  blade  rotational  velocity.  Thus 
the  vortex  system  is  assTjmed  to  be  composed  of  two  parts .  One  part  is 
formed  of  equal  geometric -pitch  helical  filaments  extending  from  each 
blade  trailing  edge  downstream  to  infinity.  The  second  part  is  a  circular 
cylinder  of  equal-pitch  helical  vortex  filaments  extending  from  the  shroud 
trailing  edge  downstream  to  infinity  and  enclosing  the  helical  sheets  shed 
from  the  blades.  There  are  the  same  number  of  helical  lines  of  intersec¬ 
tion  of  these  two  parts  as  there  are  number  of  blades . 

2 

As  was  discussed  in  Gray's  Report  ,  the  optimum  condition  requires 
3 

Horace  Lamb,  Hydrodynamics ,  Sixth  Edition,  Chap.  VII,  Dover  Publi¬ 
cations,  New  York,  19^5 « 

2 

Gray,  op.  cit.  ante  page  5* 


7 


that  in  the  ultimate  wake  this  vortex  system  must  appear  to  move  as  a 
rigid  body;  that  is,  without  distortion.  This  requirement  yields  the 
boundary  conditions  both  on  the  motion  of  the  vortex  filaments  and  on  the 
blade  bound  vortex  strength  distribution.  For  instance,  if  no  distortion 
of  the  vortex  pattern  is  to  occur,  then  the  radial  component  of  the  in¬ 
duced  velocity  must  be  zero  for  all  points  on  the  vortex  sheet  and,  in 
addition,  the  component  of  induced  velocity  that  is  normal  to  the  helical 
filaments  and  the  radial  co-ordinate  must  be  such  that  its  axial  component 
is  Identical  for  all  points  on  the  vortex  sheets .  The  induced  velocity 
component  in  the  direction  of  the  filament  is  immaterial. 

The  relationship  between  the  induced  velocity,  the  vortex  filament 
strength,  and  the  vortex  geometry  is  given  by  the  Biot-Savart  relation 
which  states  that 


where  dV.  is  the  increment  at  an  arbitrary  point  of  the  induced 

^  velocity  that  is  associated  with  a  vortex  element  of 

length  ds '  . 

K  is  the  strength  of  the  vortex  element. 

P>  is  the  angle  between  the  normal  to  the  element  ds  ' 

and  the  displacement  vector,  P  from  ds '  to  the 
point;  measured  in  the  plane  determined  by  dl'  and 

P. 

ds '  is  the  length  of  the  elemental  vortex  filament. 

P  is  the  distance  from  element  ds '  to  the  point  under 

consideration. 

The  equations  for  the  velocity  components  associated  with  an  isolated  re¬ 
entrant  vortex  filament  in  Cartesian  co-ordinates  are  given  on  page  211  of 
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Lamb's  Hydrodynamic s  and  are  repeated  below. 
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tjI  J 


(2a) 


u 
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(2b) 
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(2c) 


where,  as  shown  in  Figure  1 

x,y,z  are  the  Cartesian  co-ordinates  defining  the  point 

at  which  the  velocity  components  are  desired.  The 
z-axis  and  the  axis  of  the  vortex  system  coincide. 
The  positive  z-direction  is  in  the  direction  of 
advance  of  the  vortex  system. 

x',y',z'  are  the  Cartesian  co-ordinates  defining  the  posi¬ 
tion  of  the  vortex  element  ds ' . 


u  .u.,u_ 


are  the  induced  velocity  components  parallel  to  the 
X-,  y-,  and  z-Axes,  respectively. 


P 


(x  -  x’)^  +  (y  -  y')^  +  (z  -z 


(5) 


The  helical  filament  may  be  more  easily  described  in  cylindrical  co¬ 
ordinates  (r,  \lr,  z)  so  that  the  following  transformation  is  made. 


X  =  r  cos  xjf 
y  =  r  sin 
z  =  z 
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Coordinate  Systems . 


X '  =  r '  cos  \lf 
y '  =  r '  sin  \lf ' 
z '  =  '  +  r '  >|/ '  tan  cp ' 

s'  =  r '  itf '  sec  (p ' 


where  qj'  is  the  helix  pitch  angle  of  the  vortex  filament. 

z  '  is  the  shortest  axial  distance  between  the  x-axis  and 
the  intersection  of  the  vortex  filament  with  the  xz- 
plane . 

The  boundary  conditions  are  more  conveniently  expressed  in  temns  of  the 
velocity  components  along  and  perpendicular  to  the  vortex  sheet  filaments. 
Thus , 
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where  u 

r 


is  the  radial  induced  velocity  component . 


u 


i 


is  the  induced  velocity  component  in  the  direction  of 
the  vortex  filament  and  normal  to  the  radial  component. 


u 


is  the  induced  velocity  component  that  is  normal  to 
both  the  vortex  filament  and  the  radial  component . 


cp  is  the  helix  pitch  angle  of  the  sheet  at  the  point 
where  the  induced  velocity  components  are  to  be  com¬ 
puted  . 

Performing  the  above  transformations  yields  the  following  expressions  for 


the  Induced  velocity  components  associated  with  an  Infinite  length,  heli¬ 


cal  vortex  filament. 
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u  = 
r 


r'  z  -  z^'  -  r'  ijf'  tan  q)'  cos  (i|r'  -  it^)  J-  (^a) 


r'(z  -  z^'  -  r'  \lf'  tan  cp')  sin  (\|r'  -  \|a)  cos  (p  + 


r'  -  rr'  cos  (>!;'  _  Mf) 


.  'f}  ff 


K  r  °°  f 

"c  ‘^LiL 


r '  -  rr '  cos  (f  -  ■*1^) 


cos  q)  - 


r'  tan  cp'  [r  -  r'  cos  (t' -  t)  ]  + 

r'(z  -  z  '  -  r'  \(f'  tancp')  sin  (ijf'  -  sin  cpr  (4c) 

°  .  J  p-" 

1 

r  1  ^ 

where  P  =  ir  +  r'  -  2rr'  cos  (>]/'  -  \|f)  +  (z  -  z^'  -  r'  tan  cp’)  J  (5) 


In  order  to  obtain  the  total  velocity  components  for  the  entire 
pattern,  Equation  (4)  must  he  integrated  over  the  entire  vortex  sheet. 
As  far  as  is  known,  a  solution  to  these  equations  in  terms  of  the  ele- 


of  equal -width  segments  and  each  segment  is  replaced  by  a  finite  strength 
vortex  filament.  The  boundary  conditions  previously  prescribed  for  the 
motion  of  the  sheets  must  be  applied  to  a  finite  nvimber  of  points  on  the 
sheet  other  than  those  which  would  lie  on  the  filaments  in  order  that 
Equation  (5)  will  remain  greater  than  zero  and  the  integrands  of  Equation 
(4)  will  remain  finite.  For  this  reason,  the  calculating  points  for  this 
analysis  are  chosen  on  the  sheet  and  midway  between  the  vortex  filaments; 
i.e.,  the  points  separating  two  consecutive  vortex  segments.  As  pointed 
out  before,  the  strengths  of  these  vortex  filaments  are  unknown.  There¬ 
fore,  for  each  calculating  point,  an  equation  for  each  induced  velocltji 
component  is  obtained  in  terms  of  the  unknown  vortex  strengths.  The 
boundary  conditions  on  the  radial  and  normal  components  of  induced  velo¬ 
city  for  points  on  the  vortex  sheet  are 


u  =  0 
r 


u^  =  w  cos  cp 


(6) 

(7) 


where  w  is  the  parameter  describing  the  apparent  axial  motion  of  the 
wake  vortex  system.  For  "M"  vortex  filaments,  the  nondimens ional  equa¬ 
tions  are 
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where 
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Inspection  of  the  integrand  of  Equation  (ll)  shows  that  for  a 
boundary  sheet  strength  distribution  symmetrical  about  its  lines  of  in¬ 
tersection  with  the  inner  helical  sheets,  the  radial  velocity  due  to  the 
entire  system  is  zero  for  all  points  on  the  inner  sheets.  Therefore  the 
radial  velocity  boundary  condition  need  be  satisfied  only  for  points  on 
the  boundary  sheet.  In  addition,  for  the  pairs  of  calculating  points  on 
the  boundary  sheet  symmetrically  situated  with  respect  to  the  line  of 
intersection,  the  equations  obtained  from  (ll)  are  identical.  Thus,  the 
number  of  points  at  which  the  radial  boundary  condition  must  be  satisfied 
is  reduced  to  half  the  number  of  points  chosen  on  the  wake  boundary  vortex 
sheet  if  these  points  are  symmetrically  disposed  about  a  line  of  inter¬ 
section  . 

It  is  only  for  the  lightly  loaded  case  (i.e.,  vanishingly  small  in¬ 
duced  velocities)  that  the  assumption  of  equal  geometric  pitch  for  all 
vortex  filaments  holds.  For  cases  of  finite  loadings,  such  an  assumption 
results  in  an  incompatibility  between  the  boundary  vortex  sheet  strength 


and  the  discontinuity  in  the  normal  induced  velocity  component  u^ 
across  the  wake  boundary  sheet.  Therefore,  for  the  lightly  loaded  case. 


Equation  (l2)  need  be  satisfied  only  for  those  points  lying  on  the  inner 


helical  sheet  as  this  is  the  condition  that  determines  the  optimum  load- 


lb 


ing.  The  axial  induced  velocity  boundary  condition  should  not  be  applied 
to  the  vortex  filaments  in  the  wake  boundary  sheet.  Thus,  the  set  of 
equations  resulting  from  (12)  is  reduced  to  a  number  equal  to  the  num¬ 
ber  of  calculating  points  on  the  inner  helical  sheet. 

The  integrals  of  Equations  (ll)  and  (l2)  are  evaluated  nxmierically 
on  a  digital  computer,  yielding  a  set  of  simultaneous  equations,  one  less 
than  the  niomber  of  unknowns .  The  last  required  equation  is  discussed  at 
the  end  of  this  section. 

It  is  now  in  order  to  discuss  the  boundary  conditions  on  the  blade 
bound  vortex  strength  distribution.  Considering  first  the  inner  helical 
vortex  sheet,  the  local  strength  of  this  sheet  in  the  ultimate  wake  is 
related  to  the  change  in  blade  bound  vortex  strength  at  the  same  non- 
dimensional  radial  station  by  the  following  relation 

7  =  (f-)  (15) 

where  P  =  r(r)  is  the  blade  bound  vortex  strength  distribution. 

7  =  /(r)  is  the  vortex  sheet  strength. 

For  the  finite  strength  filaments 


“m  =  r  (r^)  -  r  (r^)  .  (fX)  dr  (14) 

“•l 

where  the  helical  sheet  between  the  radii  and  r^  has  been  replaced 

by  a  single  vortex  filament  of  strength  k  and  located  midway  between 
r^  and  r^ .  The  blade  bound  vortex  strength  must  be  zero  at  the  blade 
root.  At  the  intersection  of  the  inner  helical  sheet  with  the  boundary 


IT 


sheet,  there  can  be  no  radial  velocities  otherwise  there  would  be  a  con¬ 
tinual  distortion  of  the  vortex  pattern.  Since  the  radial  velocities  are 
zero,  there  is  no  discontinuity  in  radial  velocity  at  the  tip  of  the  heli¬ 
cal  sheet  as  the  sheet  is  crossed  so  that  by  the  Helmholtz  theorems,  the 

sheet  strength  at  this  point  is  zero.  Thus,  the  blade  bound  vortex 

cl  r 

strength  at  the  blade  tip  is  constant  or  (^~^)  =  0- 

R 

These  boundary  conditions  may  be  easily  applied  by  assuming  a 
Fourier  series  representation  of  the  blade  bound  vortex  strength  in  the 
following  nondimens ional  form. 
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(15) 


For  r  =  0:  P  =  0  so  that  A  =A,  =A^= - =  A  =  0. 

^  o  1  2  n 

—  d  r 

For  r  =  1.0:  — =  =  0  so  that  =  B,  =  B^  - - =  0. 
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Thus,  the  bound  vortex  strength  may  be  represented  by  the  infinite  series 


27TR  w  tan  cp 


■o  •  TT"  —  , 
=  B^  sin  2  r  + 


„  .  57r  -  _  .  57r  - 

B,  sin  r  +  B^  sin  r 
5  t  5  2 


r  + - 


(16) 


The  nondimens ional  vortex  filament  strength  which  replaces  the  m  seg¬ 
ment  of  the  sheet  is  then  the  difference  in  the  values  obtained  from  Equa¬ 
tion  (l6)  at  each  end  of  the  segment,  (it  is  to  be  noted  that  in  using 

Equation  (l6)  in  the  ultimate  wake,  r  =  ^  ;  while  at  the  propeller 

o 

blade ,  r  =  ^  .  ) 

I\ 

For  the  lightly  loaded  case,  there  are  no  requirements  on  the  vor¬ 
tex  strength  distribution  of  the  boundary  sheet.  It  is  to  be  noted,  how¬ 
ever,  that  the  boundary  sheet  strength  is  periodic,  having  a  period  of 
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27rR  tan  cp 
o 

- - -  .  In  addition,  it  is  symmetrical  about  its  line  of  inter¬ 
section  with  the  inner  helical  sheet  and  also  about  the  midpoint  between 
adjacent  inner  sheets.  In  a  manner  similar  to  that  of  the  inner  sheet  and 
for  computational  purposes,  one  period  of  the  boundary  sheet  is  divided 
into  equal-width  segments  and  each  segment  is  replaced  by  a  helical  vortex 
filament  of  finite  strength  and  located  at  its  midpoint. 

There  is  one  additional  important  requirement  and  that  is  that  the 
s\immation  of  the  vortex  filament  strengths  have  to  be  zero  according  to 
the  Helmholtz  theorems,  or,  in  other  words,  the  line  integral  of  the 
velocity  along  a  closed  path  completely  encircling  the  wake  must  be  zero. 
This  yields  the  one  additional  equation  required  for  the  simultaneous 


solution. 


COMPUTING  PROCEDURES 


For  the  particular  problem  under  consideration,  it  was  convenient 
to  divide  the  vortex  system  into  two  parts  and  to  further  subdivide  each 
part  into  an  equal  nmber  of  equal -width  segments .  The  inner  helical  vor¬ 
tex  sheet  was  thus  divided  into  10  equal -width  segments,  starting  at  r'  = 
0.0.  Each  segment  was  replaced  at  its  midpoint  by  a  helical  vortex  fila¬ 
ment  of  finite,  but  unknown  strength,  starting  at  r'  =0.05  and  ending 
at  the  r'  =  0.95  station  with  z'  =  0.0.  The  calculating  points  were 
chosen  at  the  ends  of  the  segments,  beginning  at  r  =  0.0  and  ending 
with  r  =  0.9-  A  portion  of  the  boundary  sheet  was  also  divided  into  10 
equal -width  segments  and  each  segment  was  replaced  by  a  helical  vortex 
filament  as  before.  Since  the  boundary  sheet  strength  distribution  will 
be  periodic  in  the  axial  co-ordinate,  the  portion  that  was  considered  ex¬ 
tended  from  the  midpoint  between  adjacent  inner  helical  sheets  and  on  the 
wake  boundary  to  the  next  midpoint  in  the  negative  axial  direction.  The 
inner  sheet  spacing  is  given  by 
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(17) 


or  nondimens ionally 

T  2ir  tan  cp 
L _ 

R  “  b 
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(18) 
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where  L  is  the  axial  spacing  between  successive  inner  helical 
vortex  sheets . 

H  is  the  geometric  pitch  of  the  vortex  filaments  forming 
the  sheet. 

b  is  the  number  of  blades  in  the  propeller. 

Thus ,  the  wake  boundary  filaments  were  chosen  to  lie  at  z '  =  ±  j 


^  0.15H  . . .  +  0.45H 


with  r'  =  1.0.  The  calculating  points  were 


0  T  TT  0  2H 

chosen  midway  between  the  filaments  as  follows:  z  =  ±  — ^ - ;  ±  — ^ —  ; 

- j  +  with  r  =  1.0.  For  the  two-bladed  propeller  wake; 

^  =  0  ;  TT  radians.  For  the  four-bladed  propeller  wake;  =  0  ;  ^  ;  tt  ; 

37r 

radians.  In  the  example  considered  in  this  report,  tan  cp^  =  I.356. 
For  given  r,  z,  and  the  integrals  of  Equations  (ll)  and 
(12)  may  be  evaluated  n\xmerically  on  a  digital  computer  for  each  vortex 
filament,  thereby  yielding  for  this  case  of  20  filaments,  kO  equations 
with  the  nondimens ional  filament  strengths  being  the  20  unknowns.  The 
intervals  A\lr"  used  in  the  Simpson's  rule  integration  were  determined 
in  the  following  manner:  First,  an  interval  of  0  to  tt  radians  was 
chosen;  the  integrands  were  evaluated  every  15°  of  this  interval,  the 
Simpson's  rule  integration  was  performed,  and  the  result  was  recorded. 

The  15°  interval  was  then  halved  to  7-5°^  "the  above  procedure  was  re¬ 
peated  and  the  result  was  compared  with  the  previous  result.  This  pro¬ 
cess  of  halving  the  interval  was  repeated  until  the  results  compared 
within  three  units  in  the  sixth  significant  figure.  When  this  occurred, 
the  next  interval,  tt  to  27r  radians  was  chosen  and  the  process  repeated. 
The  upper  limit  on  was  determined  by  comparing  the  result  of  the 
nimierical  integration  after  each  tt  interval  with  the  previous  result. 
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When  these  computations  compared  within  three  units  in  the  fifth  signifi¬ 
cant  figure,  the  integration  for  that  filament  was  terminated. 

As  was  pointed  out  in  the  previous  section,  the  calculation  of  all 
the  coefficients  appearing  in  the  sets  of  equations  was  not  necessary. 

For  the  cases  considered  herein.  Equation  (ll)  yielded  four  equations  for 
satisfying  the  radial  induced  velocity  component  on  the  wake  boundary 
(note  that  point  z  =  ±  and  r  =  1.0  was  a  point  of  helical  sym¬ 

metry  and  the  radial  velocity  was  automatically  zero),  and  Equation  (l2) 
yielded  10  equations  for  satisfying  the  normal  induced  velocity  compo¬ 
nent  on  the  inner  helical  vortex  sheet.  Due  to  the  helical  symmetry  there 
were  certain  repeating  coefficients  and  the  total  number  of  integrations 
required  were  195*  In  addition,  the  symmetry  occurring  in  the  wake 
boundary  sheet  strength  about  the  line  of  intersection  reduced  the  number 
of  unknowns  to  15-  Thus,  l4  equations  were  obtained  in  the  I5  unknowns. 
The  15th  equation  was  that  the  sum  of  the  strengths  of  the  20  vortex  fila¬ 
ments  be  zero. 

The  nondimens ional  blade  bound  vortex  strength  was  assumed  to  be 
adequately  represented  by  the  first  10  terms  of  Equation  (16).  Thus, 


27rR  w  tan  cp 
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(19) 


The  filament  strengths  were  then  simply  the  difference  in  this  equation 
evaluated  at  the  ends  of  the  vortex  sheet  segment  that  the  filament  re¬ 
placed  . 
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COMPUTER  RESULTS 


The  computer  results  are  presented  in  Table  1  and  in  Figures  2  and 

5.  In  Table  1  is  listed  the  Fourier  series  coefficients  of  Equation  (19) 

that  were  obtained  from  the  simultaneous  solution  of  the  vortex  filament 

equations  for  a  two-  and  a  four-bladed,  lightly  loaded  shrouded  propeller 

operating  at  an  advance  ratio  of  1.356.  In  Figures  2  and  3>  a-  plot  of 

Equation  (l9)  is  compared  with  the  results  of  the  potential  tank  analysis 

2 

presented  in  Gray’s  report.  The  function  K(x)  is  used  in  this 
case  for  convenience ^  where 

K(x)  =  ^  -  (20) 

'  27fR  w  tan  cp  '  ' 

o 


2 


Gray,  op.  cit.  ante  page 


5. 
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TABLE  1 


COEFFICIENTS  APPEARING  IN  THE  FOURIER  SERIES  REPRESENTATION 
OF  THE  BLADE  BOUND  VORTEX  STRENGTH  FOR  A  TWO-  AND  A  FOUR -BEADED 

SHROUDED  PROPELLER 

tan  cpQ  =  1.556 


r 

TT  —  .  5'Tr  —  . 

+  B  E  In 

27rR  w  tan  cp 

0  ^0 

—  .D-  0±I1  ^  X  T  OXli  ^  X  T  —  —  • 

X  ^  d 

+  "^19  ^  2  ^ 

Coefficient 

Two  Blades 

Four  Blades 

+  0.09089 

+  0.05305 

^3 

-  0.00297 

0.00697 

"5 

+  0.00103 

0.00099 

-  0.00048 

0.00086 

"9 

+  0.00026 

0.00028 

^11 

-  0.00014 

0.00024 

B13 

+  0.00007 

0 . 00012 

=15 

-  0.00005 

0.00007 

®1T 

+  0.00001 

0.00005 

^19 

0 

0 
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NONDIMENSIONAL  BLADE  STATION  7 


Figure  2.  A  Comparison  of  the  Results  of  Two  Methods  of  Determining  the 
Variation  in  the  Nondimens ional  Blade  Bound  Vortex  Strength 
Distribution  Function  with  Nondimens ional  Blade  Radius  for  a 
Lightly  Loaded,  Single-Rotation,  Two-Bladed  Shrouded  Propeller. 
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3.  A  Comparison  of  the  Results  of  Two  Methods  of  Determining  the 
Variation  in  the  Kondimensional  Blade  Bound  Vortex  Strength 
Distrlhution  Function  with  Nondlmensional  Blade  Radius  for  a 
Lightly  Loaded,  Single-Rotation,  Four  Bladed  Shrouded  Propeller 


The  accuracy  of  the  results  is  not  known  since  it  is  highly  depend¬ 
ent  on  the  number  of  finite  vortex  filaments  used  to  approximate  the  vor¬ 
tex  sheets .  An  estimate  of  the  accuracy  could  be  obtained  by  using  twice 


the  number  of  filaments  in  the  analysis^  but  this  was  not  considered  in 
view  of  the  extra  time  and  funds  required  for  its  accomplishment.  How¬ 
ever,  in  view  of  the  close  agreement  between  the  two  methods,  it  is  be¬ 
lieved  that  the  results  are  sufficiently  accurate  for  design  purposes. 


EVALUATION 


The  initial  comparison  of  the  blade  bound  vortex  strength  distribu¬ 
tion  for  the  two-bladed  shrouded  propeller  as  determined  by  the  computer 
analysis  with  that  measured  on  the  potential  tank  model  showed  a  differ¬ 
ence  of  about  seven  per  cent  at  all  blade  stations,  the  computer  quanti¬ 
ties  being  the  higher.  This  was  believed  to  be  excessive  in  view  of  the 
estimated  accuracy  of  the  analyses  so  that  a  re-evaluation  of  the  poten¬ 
tial  tank  data  was  undertaken.  A  small  computational  error  was  discov¬ 
ered  in  the  quantity  that  was  used  to  nondimensionalize  the  potential  tank 
measurements  for  the  two-bladed  propeller  case.  These  results  were  cor¬ 
rected  and  are  plotted  in  Figure  2  for  the  comparison  purposes.  A  similar 
check  of  the  four-bladed  propeller  case  did  not  show  a  similar  error,  so 
that  this  data  is  as  presented  before. 

The  comparison  between  the  computed  nondimens ional  blade  bound 
vortex  distribution  and  the  measurements  obtained  from  the  potential  tank 
model  is  very  good  for  the  two-bladed  shrouded  propeller.  As  may  be  seen 
from  Figure  2,  the  difference  is  about  2  per  cent  below  the  mean  mea¬ 
sured  value  at  the  blade  tip.  Since  the  potential  measurements  were  made 
at  the  center  of  a  two -turn  helical  model,  it  was  felt  that  the  end 
effects  of  the  model  might  be  causing  this  small  difference.  Accordingly, 
the  integrations  in  the  computer  method  were  limited  to  the  same  length 
and  the  Fourier  series  coefficients  were  re -determined.  It  was  found  that 
the  effect  was  in  the  right  direction  but  was  very  small,  bringing  about 
an  Increase  of  0.2  per  cent  in  and  about  0.5  pei"  cent  in  .  There- 
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fore,  Inasmuch  as  the  estimated  measurement  error  for  the  potential  tank 
results  was  ±  2  per  cent  of  the  maximum  readings,  it  was  felt  that  the 
agreement  between  the  two  methods  was  sufficient  and  further  inquiry  into 
this  small  difference  was  not  warranted.  However,  the  difference  may  be 
attributed  to  small  errors  in  model  construction  and/or  to  the  approxima¬ 
tion  introduced  into  the  computer  method  by  replacing  a  vortex  sheet  with 
a  finite  number  of  finite  strength  vortex  filaments . 

As  may  be  seen  from  Figure  3^  the  difference  in  the  computed  blade 
bound  vortex  distribution  and  the  measurements  from  the  potential  tank 
model  for  the  four-bladed  shrouded  propeller  is  also  about  2  per  cent  at 
the  blade  tip,  but  is  considerably  more  near  the  blade  root.  Inasmuch  as 
the  four-bladed  potential  model  was  much  more  difficult  to  construct  and 
therefore  of  less  accuracy  than  the  two-bladed  model,  it  is  believed  that 
this  increased  percentage  difference  is  due  primarily  to  model  construc¬ 
tion  errors  and  to  the  small  warpage  that  occurred  in  the  model  as  it 
absorbed  the  electrolyte . 

In  view  of  the  success  achieved  with  the  computer  method  and  from 
past  experience  of  the  difficulties  encountered  in  the  construction  of 
potential  tank  models  of  the  wakes  of  shrouded  propellers,  it  is  con¬ 
cluded  that  the  computer  method  is  the  more  efficient  means  of  obtaining 
the  desired  results. 

From  Table  1,  it  may  be  seen  that  for  the  cases  investigated,  the 
first  four  coefficients  of  the  series  representation  of  the  blade  bound 


vortex  distribution  are  sufficient  for  engineering  design  purposes.  The 
use  of  this  form  of  representation  has  two  primary  advantages.  The  first 
has  been  discussed  previously  and  allows  the  boundary  conditions  on  the 
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diFtrilDution  to  be  expressed  more  easily.  The  second  advantage  would 
appear  in  the  presentation  of  the  results  if  the  present  work  were  to  be 
extended  to  generate  the  necessary  optimiim  parameters  for  design  purposes 
by  performing  the  computations  for  a  range  of  advance  ratios  and  numbers 

4  1 

of  blades.  Previous  analyses  for  free  and  ducted  propellers  have  pre¬ 
sented  a  family  of  curves  or  a  series  of  tables  for  the  nondimens ional 
blade  bound  vortex  distribution  versus  nondimensional  blade  radius  (simi¬ 
lar  to  Figures  2  and  3)  for  various  values  of  a  parameter  \  which  is  the 
tangent  of  the  helix  angle  of  the  blade  tip  vortex  filament .  For  values 
of  K  other  than  those  presented  in  the  reports,  an  interpolation  must  be 
performed.  However,  if  the  sine  series  representation  is  used,  it  is  be¬ 
lieved  that  the  coefficients  may  be  plotted  as  a  function  of  advance  ratio 
and  an  explicit  expression  of  four  or  five  terms  may  then  be  obtained 
which  will  yield  the  required  distribution  versus  blade  station  by  a  rela¬ 
tively  simple  calculation.  This  will  reduce  the  number  of  curves  or 
charts  required  for  adequate  data  presentation  and  should  also  require  a 
fewer  number  of  advance  ratios  for  which  the  complete  computer  analysis 
must  be  run.  This  is  a  particular  advantage  in  that  about  8  to  10  hours 
of  computing  time  is  required  for  each  run.  It  thus  appears  that  a  very- 
real  savings  could  be  obtained  by  use  of  the  series  representation. 

The  analysis  presented  herein  is  applicable  only  to  lightly  loaded 

shrouded  propellers.  In  the  case  of  the  free  propeller,  the  lightly 

4 

loaded  analysis  was  extended  to  the  heavily  loaded  case  by  a  simple 
4 

T.  Theodorsen,  Theory  of  Propellers,  McGraw-Hill  Book  Company, 
Inc.,  New  York,  19^8. 

^achmindji,  op.cit.  ante  page  5* 
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redefinition  of  one  of  the  parameters.  The  extension  to  the  heavily  loaded 
shrouded  propeller  does  not  appear  to  he  as  simple  because  of  the  addi¬ 
tional  compatibility  requirement  between  the  inner  helical  vortex  sheet, 
the  boundary  vortex  sheet,  and  the  discontinuity  in  induced  velocity  as 
the  boundary  sheet  is  crossed.  This  requirement  was  neglected  in  the 
present  analysis  and  therefore  restricts  the  results  obtained  from  this 
and  similarly  restricted  analyses  to  those  cases  in  which  the  induced 
velocity  is  small  in  comparison  with  the  free-stream  velocity.  If  there 
is  no  urgent  need  for  the  lightly  loaded  design  parameters  in  the  near 
future,  it  is  suggested  that  this  approach  be  discontinued  and  any  further 
efforts  be  directed  toward  an  analysis  of  the  heavily  loaded  case,  using 
the  same  methods  presented  herein  but  including  the  additional  condition 
on  the  wake  boundaiy  vortex  sheet. 

It  was  felt  that  a  comparison  between  the  method  of  this  report  and 
that  of  Tachmlndji^  would  be  worthwhile.  Therefore,  the  nondimens ional 
blade  bound  vortex  strength  distribution  for  a  four-bladed  shrouded  pro¬ 
peller  operating  at  an  advance  ratio  of  one-third  was  calculated  and  the 
comparison  is  shown  in  Figure  4.  It  was  expected  that  better  agreement 
between  the  two  methods  would  be  obtained.  Limitations  on  time  and  funds 
did  not  permit  further  investigation  of  the  differences  shown. 

As  an  additional  comparison,  the  distribution  for  an  Infinitely 
bladed  shrouded  propeller  was  calcula,ted  using  an  expression  given  by 
Theodorsen.^  This  is  also  shown  in  Figure  4. 


^achmind  j  i , 


Theodorsen, 


op .  c it .  ante  page  5 • 
op.  cit.  ante  page  JO* 
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